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SECTION I 
INTRODUCTION 


This report describes the aerodynamic design of a tran- 
sonic axial-flow compressor inlet stage designed for a high 
flow rate per unit frontal area and relatively high aerodynamic 
loading. The performance objectives of this stage were 
derived from a preliminary design study of a multi-stage 
compressor for an advanced turbojet engine. In the course of 
this study it became apparent that the mast serious aerodynamic 
design problems of the overall compressor were associated 
with the first stage. Specifically, it was necessary to 


‘design the first stage with relatively high diffusion factors 


and stator hub Mach number in addition to high flow per 

unit frontal area in order to keep loading levels in the 
remaining stages within reason while maintaining a combustor 
inlet Mach number less than 0.4. The design criteria finally 
chosen for the first stage were felt to represent the best 
compromise from the standpoint of the overall multi-stage 
compressor. This stage was chosen as the object of an 
independent research program because it presses the present 
state-of-the-art with respect to nearly all of its design 
parameters, Also, it appeared to be a suitable test vehicle 
on which to assess the usefulness of several vortex generator 
configurations which have recently proven successful ona 
more highly loaded supersonic compressor stage. 


The preliminary design is discussed in the second 
section of this report. The third section fully describes 
the detailed aerodynamic design assumptions and procedures. 
The final aeredynamic design results are. presented in the 
fourth section. This includes a complete aerodynamic 
description of the compressor stage in addition to 
specification of the guometry. 
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SECTION II 
PRELIMINARY DESIGN 


1. CRITERIA 


All of the criteria defining the basic parameters of this 
compressor stage resulted from the design goals of a hypothet- 
ical turbojet engine of which this compressor comprised the 
first stage. A flow per unit frontal area of 39.7 lb/sec per 
square foot was established at the outset as a design goal. 

A corrected tip speed for the first stage of 1500 ft/sec was 
also established early in the design study as the maximum 
value consistent with turbine stress considerations for this 
particular application. Most of the rest of the compressor 
characteristics resulted from the overall objective of 
designing a compressor for a turbojet engine with the minimum 
number of stages consistent with an overall isentropic effi- 
ciency of 0.84, a specified pressure ratio, and a combustor 
inlet Mach number less than 9.4. Further details concerning 
design criteria for the complete multi-stage compressor have 
been deliberately omitted from this report to avoid restric- 
tions imposed by security classification, 


2. PROCEDURE 


The preliminary design of the multi-stage compressor was 
accomplished with the computer program described in 
Reference 1, This computer program performs an axisymmetric, 
full radial equilibrium analysis of the compressor flow 
field using the streamline curvature solution technique. 
Through a series of iterations, it attempts to maximize the 
performance of each successive stage by driving the design 
teward one or more of a number of specified limits and with a 
number of specified constraints. The racial distribution 
of energy addition in each rotor is specified by the user 
as a non-dimensional total pressure distribution, The aerody-~ 


namic parameters for which limits must be supplied for each 


stage are: — 


(1) Rotor tip diffusion factor 
(2) Stater hub diffusion factor 
(3) Stator hub Mach number 

(4) Rotor hub relative exit angle 
(8) Rotor tip exit whirl velocity 


The first three of these limits were found most useful for 
this investigation. The fourth and fifth were simply set at 
values which would avoid constraining the results. The general 
precedure followed was te start by specifying relatively 
conservative values for these limits This initially 
resulted in a design having insufficient pressure ratio and 
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higher efficiency than required within the maximum number of 
stages allowed. Subsequently, the controlling limits were 
gradually raised and some adjustments were made to the radial 
distributions of rotor work and mid-streamline axial velocity 
ratios until the desired overall pressure ratio was achieved at 
the objective efficiency and with what appeared to be a good 
balance of conditions throughout the compressor. In the 
first part of the preliminary design, the program option 
involving specified mid-streamline axial velocity ratios was 
employed and the annulus walls were allowed to float within 
limits. Final fine tuning of the preliminary design was 
accomplished with the annulus geometry frozen. A minor modi- 
fication was made to the computer program for convenience in 
arriving at this design. This is described in Appendix A. 


3. LOSS ASSUMPTIONS 


The losses attributed to each blade element and expressed 
as a relative total pressure loss coefficient have been 
assumed to be equal to the sum of two components; one 
associated with diffusion occurring in the profile houndary 


‘layers and one related to the presence of shock waves in each 


blade passage. Shock losses were estimated according to the 
familiar Miller-Lewis-Hartmann model described in Reference 2. 
This consisted of assuming the shock-related portion of the 
relative total pressure loss coefficient to be that 

resulting from a normal shock, the upstream Mach number of 
which was equal to the arithmetic average of the relative 

inlet Mach number and the suction surface Mach number at 

the shock impingement point. This latter Mach number was 
assumed te differ from the rolative inlet Mach number by a 
Prandtl-Meyer expansion through a specified number of degrees 
of turning. The previously mentioned modification (presented 
in Appendix I) to the computer program described in Reference 1 
consisted of an option allowing the program user to specify 

the supersonic turning angle. This proved to be 

more convenient for this design exercise than other options 
available to accomplish the same purpose. The final version 

of the preliminary design assumed the supersonic turning angle 
to vary linearly from 12 degrees at the hub to 2 degrees at 

the tip for the reter and assumed a constant 15 degrees for the 
stator. These assumptions were varied somewhat in the aarly © 
stages of the preliminary design until the approximate blade 
geometry became evident. Where the Mach number relative to a 
blade element is less than one, the procedure was handled 
Slightly differently in order to effect a smooth transition in 
the loss distribution. The suction surface Mach number was 
assumed to be that resulting from a Prandtl-NMeyer expansion from 
Mach 1 through the given supersonic turning angle. Then, the 
Upstream shock Mach number — assumed equal to 


My = ele (1 + Mss) 


instead of the previously used average value. Whenever Mx was 
less than or equal to 1.0, shock loss was assumed to be zero. 


Losses resulting from diffusion were predicted in the 
Lieblein manner by a relative total pressure loss parameter 
versus diffusion factor relatiorship. The computer 
program of Reference 1 allows this relationship to be defined 
at 10, 50 and 90 percent span, independently for each blade 
row. lnitially curves taken from Figures 33 and 34 of 
Reference 3 were employed as representing the latest published 
revision of the original correlation. However, as the pre- 
liminary design began to take shape, the objective pressure 
ratio was achieved within the desired number of stages, but 
with efficiencies bordering the incredible, several points 
above the best ever obtained under similar or easier circum~ 
stances. Even greatly increasing the supersonic expansion 
angles to clearly impossible values did not bring the predicted 
efficiencies within the realm of credibility. A closer 
examination of the loss parameter correlation of Reference ? 
led to the conclusion that the prediction shown for rotor 
profile loss appeared extremely optimistic for diffusion 
factors less than 0.5. representing essentially the entire 
range of interest for this design. The floor level of the 
loss parameter curves had been shown as approximately 0.002. 
For this design, this was raised to 0.009 and new curves were 
created which passed through the values given by Reference 3 
at a diffusion factor of 0.7, and were extrapolated to 1.0 
for the rotor, and retaining the prediction of Reference 3 for 
the stator, reasonable results were achieved. The distributions 
used are reproduced in Figures la and lb of this report. 


4, DESIGN PHILOSOPHY 


a. Velocity Triangles 


The choice of velocity triangles throughout the multi- 
stage compressor was moet strongly influenced by the desire to 
achieve the maximum pressure ratio per stage consistent with 
reasonable off-design performance and the design efficiency 
. obtective. Since all rotors were expected to operate at 
relatively high levels of diffusion, high rotor relative Mach 
numbers and the associated shock losses were viewed as an 
acceptable penalty for the resulting performance. Although 
vaplable inlet guide vanes (IGV) might be desirable to 
expand the part-speed operating envelope, they were not 
desired from the point of view of altering the relative inlet 
Mach number of the first rotor, nor could significantly 
. cambered IGV be tolerated at design point operation because 
of the high flow rate and r-sultant danger of choking. 
Consequently, a design approach was adopted which employed no 
IGV, assumed that the first stage stator would be a rigid 
structure supporting the front bearing and transmitting services, 
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and assumed that satisfactory off-design performance could be 
obtained by variable stators in the second and later stages, 
and an articulated trailing edge in the first stator and/or 
bleeds if necessary. Furthermore, all stators were designed 
to turn the flow back to the axial direction. 


b. Axial Velocity Ratio 


Three factors guided the magnitude and distribution of 
axial velocity ratios selected for each blade row. First, the 
entire multi-stage compressor required an overall axial velocity 
ratio less than unity in order that the combustor inlet Mach 
number not exceed u.40, one of the initial conditions. Rather 
than force just one or two stages to accept a particularly low 
axial velocity ratio, it was decided to allow the stage 
axial velocity ratios co gradually decline over the first 
few stages and then remain approxinately constant over the 
balance of stages. A mid-streamlire axial velocity ratio of 
approximately 0.98 was chosen for the first stage. The 
second two factors cetermined the distribution of axial 
velocity ratio between rotor and stator. Rotors have generally 
proven capabl of operating at nigher loading levels than 
stators havir~ «..ial losses. This favored axial velocity 
ratios less tnau unity in rotors and greater than unity in 
stators. Also, as the axial velocity ratio across a rotor is 
reduced, the camber required to produce a given change in 
swirl velocity is also reduced. Consequently, since average 
rotor relative Mach numbers were much higher than average 
stator Mach numbers, low rotor axial velocity ratios tended to 
minimize supersonic turning angles, and consequently shock 
losses for a given level of solidity. The mid-streamline 
values ultimately arrived at for the first stage were 0.815 
for the rotor and 1.206 for the stator. 


ec. Annulus Shape 


The shape chosen for the compressor annulus represented 
a compromise between maintaining high enough blade speed to 
facilitate obtaining high pressure ratio per stage while at the 
Same time not permitting the hub/tip radius ratio at the 
compressor exit to become too high. There were also mechan- 
ical considerations involved such as ease of maintaining tip 
clearance and stresses in the last rotor. The shape. 
ultimately chosen incorporated a constant outer diameter over 
the first three blade rows. The oviter diameter was then 
reduced somewhat across each stator starting with the second 
Stage. All rotors downstream of the second maintained a 
cylindrical outer diameter for ease in controlling tip 
clearance, although these diameters were progressively 
reduced, This choice of contour produced an adverse 
aerodynamic condition with respect to the effects of 
streamline curvature on aerodynamic loading. However, this 
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aerodynamic effect had such small impact on loading levels and 
losses that it was considered a desirable trade-off against 
better control of tip clearance. 


d. Aspect Ratio 


The choice of aspect ratio was governed by four factors, 
not all of which were always important in any one stage. In 
the first stage, the most important factors were hub ramp 
angle, mechanical stresses, and aerodynamic stability. High 
aspect ratio tended to minimize stage weight and rotor root 
stresses but tended toward steep hub ramp angles, greater 
flutter sensitivity necessitating part-span dampers, and less 
stall margin. Lower aspect ratio improved these last three 
conditions at some sacrifice in the first two, but if carried 
tco far, could lead to unacceptable rotor root stresses due to 
increasingly acute angles between rotor blade and hub near 
the trailing edge. Downstream of the first stage, ramp angle 
and blade stresses became progressively less important and a 
fourth factor, manufacturing cost, became more important. 

Lower cost favored lower aspect ratio in later stages. The 
computer program of Reference 1 employed an aspect ratio based 
upon the ratio of leading-edge span to the axial depth of a 
blade row. Using this definition, a value of 2.0 was selected 
for the first rotor, which resulted in an initial hub ramp 
angle of approximately 30 degrees and appeared to offer reason- 
able prospects of not requiring mid-span dampers. The aspect 
ratio of each of the next several blade rows was chosen strictly 
on the basis of adjusting the axial length of each blade 

row such that the hub ramp angle declined monotonically and 
smoothly from the inlet value to nearly zero. An aspect ratio 
of 1.° was chosen for the last few blade rows as representative 
of the minimum value offering an acceptable ratio of blade span 
to circumferential spacing. 


5. RESULTS 


The final results of the preliminary design pertaining 
to the first compressor stage and its immediate environment 
are shown in the following pages of computer printout. The 
information presented is sufficiently complete that, with the 
aid of Reference 1 and the associated computer program 
(available from COSM1C), anyone so inelined should be able to 
duplicate the results. 


One of the input requirements of this conputer program was 
a nondimensional spanwise total pressure distribution in the 
exit plane of each rotor. For the first few design iterations, 
total pressure leaving each rotor was assumed constant as a 
matter of convenience. As soon as the design began to evelve, 
however, it became necessary to increase the total pressure 
leaving the hub of each of the first few rotors in order to 
compensate for stator hub losses and to obtain a reasonable 
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balance of conditions throughout the compressor. 


; The diffusion factor at the rotor tip proved to be the 

limiting parameter in the first stage. This reached the specified 

; limit of 0.52. The resulting Mach number relative to the 

: stator hub was barely supersonic at about 1.02 and dropped 
rapidly with increasing radius. The stator hub diffusion 

: factor was somewhat high at about 0.55, but the predicted 

: iosses were acceptable and this appeared to be a necessary 

compromise for the benefit of subsequent stages downstream. 
The mid~streamline axial velocity ratio of 0.815 across the 

: first rotor kept the relative turning angle across the rotor 

: tip down to 2.4 degrees. This was thought to offer good 

: prospects of minimizing rotor shock losses and deviation 

: angles. The solidity variation in the first rotor assumed 

i for the preliminary design (a linear variation from 

: approximately 2,0 at the hub to 1.3 at the tip) proved 
eventually to be somewhat low. This solidity was increaesed 
during the detail design procedure discussed later in this 
report. 


The performance predicted for the first stage was a total 
pressure ratio of 1.93 at an isentropic efficiency of about 84 
percent. Reducing the aerodynamic loading would have led to 
sligh.ly higher efficiency and sutstantially less design 
risk in this staze. However, it would have required higher 
loading in downstream stages and/or a higher discharge 
Mach number from the lest stage, neither one of which was 
viewed as desirable fom ar. engine design standpoint. 


A matter of practical interest to users of the computer 
program associated with Reference 1 concerns the tolerances 
applied to portions of whe iteration scheme. Most of the 
tolerances recommended in Reference 1 were found satisfactory, 
However, for this denign, much less difficulty was 
experienced with the program when the loading limit tolerarce 
Was increased from the recommended 0.01 to 0.033. 
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SFCTION ITI 
DETAILED AERODYNAMIC DESIGN 


1. COMPUTATIONAL METHOD 


The detailed aerodynamic design of the first compressor 
stage was accomplished with an early version of the computer 
program described in Reference 4, and employing the "stream- 
line curvature" method of computation. Although conceptually 
similar to the method employed for the preliminary design, 
much greater precision was incorporated into the detail 
design by the addition of many more axially-distributed 
computing stations to define the flow path, including four 
stations within the interior of each blade row. The detail 
design program permitted the use of curvilinear computing 
stations, and the se were adjusted periodically to insure a 
high degree of coincidence with blade-row edges. Twenty-one 
streamlines \.:.e employed for the detail design in contrast 
to the eleven employed for the preliminary design. 

seicene iets 

A solution was obtained through an iterative numerical 
procedure according to which the equations of momentum, 
continuity, and energy are simultaneously solved at each 
computing station in sequence throughout the compressor, 
After each "pass" through the compressor, the streamline 
geometry is updated and the procedure is repeated until the 
cuénges occurring from one pass to the next fall within a 
specified tolerance, whereupon the solution is considered to 
be "converged. " The most significant assumptions are that 
the flow is axisymmetric and can be described by a series of 
concentric streamsurfaces <icre?s which no mass or 
momentum is transferred. Within this framework, the "full" 
radial equilibrium version of the momentum ecuation is 
satisfied at each streamsurface/computing-station intersection 
whereby effects of streamiine curvature and entropy gradients 
are included. For the computations made within a blade row, 
blade force terms are included in the momentum equation 
in the form of a body-force field assumed to act in a 
direction everywhere noyrmal to the three-dimensional surface 
formed by the stacked camber lines of each blade row. A 
correction for the meridtonal entropy gradient was anonudes 
according to Reference 5, 


An allowance for boundary Layer blockage was incorporated 
to account for boundary layer development on both blades and 
annulus walls. The absolute value of blockage was estimated 
according to a Simplified formula for a turbulent boundary 
layer on a flat plate, with meridional velocity and length 
substituted for absolute velocity and path length. The 
methoe used is substantially as described in Reference 6 and, 
although relatively crude, it has produced reliable results 
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uader a variety of circumstances. The resulting total 
blockage calculated at cach computing station is linearly 
distributed from hub to tip as a function of radius and is 
introduced as a factor in the continuity equation. Within 
the interior of blade rows, the additional blockage resulting 
from the finite thickness of the airfoils is also taken 

into account. 


2. ATRYOIL SELECTION AND OPTIMIZATION 


When a design is accompiished involving computing stations 
internal to a blade row, some optimization criteria mvst.be 
specified according to which the most efficient airfoil 
geometry can be selected from among a variety of possible 
designs satisfying the same end conditions. The axial 
distribution of static pressure along each streamsurface, as 
computed by the axisymmetric flow analysis, was selected as 
the most appropriate parameter to optimize for this design. 
This was chosen on the basis of being the parameter most 
closely related to blade surface boundary layer behavior 
which could be calculated with some degree of accuracy. The 
minimum staiic pressure gradient along each streamsurface 
would obviously occur with a linear variation between leading 
and trailing edge. Since this condition could probal.iy not be 
realized in practice, high deviation angles and losses would 
be the probable result. The "optimum" 4 axial distribution of 
(circumferentially averaged) static pressure along each 
streamsurface has been defined as one which is approximately 
Jinear over the first three quartezs of a blade row and then 
declines smoothly to nearly zero at the trailing edge in 
deference to tne Kutta condition. 


Two basic design approaches are feasible in conjunction with 
the optimization criteria described above. In the more 
traditional approach, one can assume the blade geometry, 
solving for the equilibrium flow field using specified 
relative flow sigles as input to the aerodynamic program. In. 
this design approach, the parameters defining the blade 
geometry are adjusted through a series of iterations until the 
optimization criteria are achieved over as much of the 
blade surface as possible. This technique has two shortcomings. — 
First, one does not have infinite flexibility to adjust the 
shape "of blades of a specifie geometric family, so the 
optimization criteria can rarely, if ever, be achieved over 
the full span of a blade. Second, the use of specified 
relative flow angles as input to the streamline«curvature- 
type aerodynamic analysis program can lead to numerical 
instabilities anc convergence difficulties in the calculation 
procedure at high subsonic through-flow Mach numbers. The 
design approach using arbitrary airfoils avoids both of these 


difficulties and was chosen for this design. Using this 


technique, the designer must assume the work distribuvion 
(total temperature or enthalpy) along streamlines through the 
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rotor and the swiri velocity (or preferably radius times 

swirl velocity) distribution along streamlines through the 
stator, The aerodynamic flow field analysis then prceduces a 
set of relative flow angles to which airfoils must be matched. 
The technique and computer program described in Reference 7 
were developed for this purpose. This procedure is iteratively 
repeated until the optimization criteria are met over the 

full span, airfoil metal blockages used in the aerodynamic 
analysis are mutually consistent with thoge calculated by 

the biade generation program, and the computing stations 

used to represent blade leading and trailing edges are a 

close match to the envelope of the stacked airfoil. In 
pursuing such a design approach, it is absolutely essential to 
insure that the distributions of work through rotors and 
radius-times-swirl-velocity through stators vary smoothly 

both along streamlines as well as along computing stations. 
Any lack of smoothness in these distributions will be 

directly translated into peculiarly shaped airfoils having 
undesirable mechanical and propably undesirable aerodynamic 
properties. A point of practical interest to designers is 
that linear distributions of these variables along streamlines 
provide an excellent starting point for a design. The 
adjustments away from the linear distribution required to 
achieve the optimization criteria are usually not very large. 


3; INTRA-BLADE AERODYNAMIC ASSUMPTIONS 


In conjunction with a through-the-blade-row design technique 
such as described above, some assumption must be made for 
three critical meridional distributions within each blade 
row. These are deviation angle, blockage, and losses, listed 
in decreasing order vf importance. At high relative flow 
angles such as encountered near the tip of a high-speed rotor, 
the performance is extremely sensitive to changes in 
relative fiow angle. There is virtually no experimental 
data available to define a reliable correlation applicable to 
the interior region of a compressor blade row. Furthermore, 
few if any reliable analytical techniques are currently 
available to predict circumferential average relative flow | 
angles through a transonic blade row. Consequently, distribu- 
tions were assumed which met the following criteria: 


(1) Deviation at the leading edge must equal 
the incidence angle — ; 


(2) Deviation at the trailing edge must 
equal the value predicted by conventional 
empirical deviation angle correlations 


(3) Deviation angles in the covered 
portion of the passage should be extremely 
small 


i (4) The rate of increase of deviation angle 
approaching the trailing edge must approxi- 
Be : mately equal the rate of change of camber in 

that region in order to approach a Kutta 
la condition at the trailing edge. 


The distributions chosen for this design are vepeneed in 
Figure 2. . 


The absolute level of blockage at each computing station 
was calculated according to the crude boundary layer model 
mentioned earlier in Section III.1. Although the 
computation was originally conceived as an annulus wall 
blockage model, it was employed here as a wake blockage 
model and linearly distributed spanwise. The ratio of 
hub blockage to tip blockage within each blade row was set 
equal to the ratio of hub solidity to tip solidity. The 
distribution is shown in Figure 3. 


The aerodynamic calculation was least sensitive to the 
meridional distribution of losses within each blade row, 
Also, viewed in the meridional plane, the rotor shock 
system covered most of the axial depth of the rotor, 
Consequently, total pressure losses were simply linearly 
distributed through each blade row. 


4, ANNULUS SHAPE 


When designing a compressor incorporating blades of a 
specified geometric family, e.g. multiple circular arc, etc., 
coupled with a through-the-blade-row design approach, it is 
usually found necessary to adjust the annulus wall contour 
within the blade passage in order to satisfy optimization 
criteria. This frequently results in locally undesirable 
pressure gradients due to streamline curvature effects near 
the walls, With arbitrary airfoils, such as employed for this 
: design, this situation can be largely avoided since the 
4 airfoil can adjust to meet optimization objectives over its 
entire surface, independently of the wall contour. Conse- 
quently, a specific objective in laying out this design was to 
maximize the radii of curvature defining the hub flowpath. | 
As a result, the spinner contour is a circular are tangent to 
the rotor hub, the rotor hub was made conical, and the hub 

a flowpath from the rotor trailing edge to the exit plane is a 

4 single circular arc tangent to the rotor hub and to a cylinder 
wf about half a stator chord length downstream of the stator 
4 exit plane, A reflex curvature in this contour was deliberately 
we 


avoided. The tip flowpath consisted of a circular-are 
bellmouth tangent to a cylinder which extended past: the praee 
exit plane. 
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5. ROTOR ASPECT RATIO 


Whereas the preliminary design employed only a single, 
radial computing station to represent both the trailing 
edge of one blade and the leading edge of the following blade, 
the detail design used individual curvalinear stations to 
represent each edge as well as intermediate locations. The 
axial depth of the rotor at the hub was selected as the more 
important axial dimension to preserve without change from the 
preliminary design, since it directly affected rotor hub 
ramp angle. During the course of detail design iterations, 
it was necessary to reduce the inlet hub/tip radius ratio 
from the preliminary design value of approximately 0.35 
to 0.31 in order to eliminate choking problems. At the 
rotor trailing edge, the hub/tip radius ratio was similarly 
reduced from 0,54 to 0.52. These hub changes caused the 
final ramp angle to rise to 32.5 degrees, the axial depth 
having remained constant. Final rotor aspect ratio was 1.32 
based upon mean rotor span and the average of the chord 
length of blade sections on streamlines 1, 11, and 21. 


6. STATOR ASPECT RATIO — 


The stator aspect ratio was decreased somewhat from the 
preliminary design, primarily as a result of the choice of 
hub flowpath contour described in Section III.4. The stator 
exit-plane area was preserved from the preliminary design, 
fixing the hub radius at this location. This radius inter- 
sected the circular arc defining this contour about twenty 
percent farther downstream from the rotor than in the prelimi- 
nary design. The resultant final stator aspect ratio was 
1.255, based upon mean blade span and the average of the chord 
length of blade sections on streamlines 1, il, and 21. 


7. ROTOR SOLIDITY AND THICKNESS 


Rotor solidity and thickness are discussed together 
because of their relationship from a structural standpoint. 
Maximum section thicknesses of 6.0 percent chord at the hub 
and 2.5 percent chord at the tip were selected as the minimum 
values likely. to be acceptable structurally. This subse- 
quently proved adequate. Preliminary centrifugal stress 
calculations also indicated that rotor airfoil cross- 
sectional area should increase by a factor of approximately 
two from tip to hub. Rotor hub chord was already established 


within narrow limits on the bases of flow area contraction 


ratio and ramp angle, as discussed in Section III.5. 

These constraints fixed the tip chord length at approximately 
4.0 inches. Aerodynamic considerations dominated the choice 
of tip solidity, with no constraint placed upon maximum hub 
solidity. It was desired that the tip solidity be high 
enough so that a weak oblique passage shock extending from 


23 


ee eee 


Gre aetenes 


aE Te Ne gine penipreeteneig OF mutate lap ate 


a ae 


the leading edge of one airfoil to the suction surface of 
another would be fully captured within the passage. This was 
felt to enhance aerodynamic stability. A tip solidity of 
approximately 1.5 was chosen so that such a shock would 
intersect the suction surface at about 90 percent chord. 

This tip solidity plus the chord length set the number of 
rotor blades at 20. An even number was chosen for convenience 
in balancing the rotor. The resultant hub solidity was 
approximately 3.2. Because of the high shock losses which 
were potentially possible near the tip and the danger of 
choking near the hub, leading edge wedge angles were minimized 
by locating the position of maximum airfoil thickness relatively 
far aft. A linear spanwise distribution was specified, 
varying from 56 percent chord at the hub to 70 percent 

chord at the tip. The airfoil leading edge radius was 
initially fixed at 0.005 inch, spanwise constant, in order to 
minimize shock losses, However, high leading edge stresses 
near the hub subsequently caused the leading edge radius to 

be flared to triple this value within the last fifty percent 
of span approaching the hub. This distribution is shown in 
Figure 4. During the first few iterations of the detail 
design, linear spanwise distributions of chord length and 
thickness were specified. However, the resultant solidity at 
about two-thirds span appeared to be dangerously low from an 
aerodynamic standpoint. The final configuration was 

achieved by specifying the solidity distribution shown in 
Figure 5, increasing the chord length and decreasing the 
thickness near two-thirds span to preserve the satisfactory, 
nearly linear area distribution. The final spanwise distribu~ 
tions of thickness-to-chord-ratio, and streamsurface-section 
and cartesian-section area are presented in Figures 6 through 8 
respectively. 


8. STATOR SOLIDITY AND LEADING-EDGE SWEEP 


The stator solidity distribution was defined purely on the 
basis of aerodynamic considerations related to loading, | 
operating range, and leading-edge sweep. Sweep was used at 
the hub to minimize shock losses and to maximize incidence | 


. range. The locally high hub solidity also reduced the 


diffusion factor at the hub relative to the preliminary 
design, which had shown an increase in that region. Chord 
length was approximately fixed according to the considerations 
related to aspect ratio and hub flowpath contour described in 
Section ItI.6. An odd number of blades having no common 
divisor with the 20 rotor blades was desired to minimize the 
chance of exciting any resonant frequencies in the rotor 
through rotor-stator interaction. A minimum solidity of 
approximately 1.6 was selected because of the high turning 
required of the stator (varying from 40 to 50 degrees) and 
the high average Mach number level varying from about 0.7 to 
1.0, These conditions were satisfied by 31 stator blades. 
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Initially, the stator trailing edge was a radial line 
coincident with the stacking axis. However, this resulted in 
a tip solidity of about 1.48 which appeared risky in relation 
to the turning, Mach number, and incidence variation which 
this region would experience over the expected operating range. 
Tip solidity was increased to 1.62 by slanting the trailing edge 
linearly aft from hub to tip with respect to the stacking 
line. The leading edge at the hub was located as far 

forward as mechanical clearances permitted, resulting in a 
hub solidity of 2.8. The angle of sweep was chosen such that 
the component of Mach number normal to the leading edge was 
approximately 0.4, the Mach number below which no further 
inerease in low-loss incidence range has generally been 
observed in cascade experiments. The sweep with respect to 
the approach flow was smoothly reduced to zero at a radius 
near mid-span where the Mach number had dropped to about 0.85. 
The leading edge then swept slightly forward toward the tip 
simply to maintain adequate solidity at the outer radii. 

The final chord length and solidity distributions are shown 
in Figures 9 and 10. Maximum thickness was located at fifty 
percent meridional chord and its absolute value varied 
approximately linearly from 4.0 percent chord at the hub to 
6.0 percent chord at the tip. Leading and trailing edge 
radii were approximately constant at about 0.0055 inch. 


9. INCIDENCE ANGLE 


The use of arbitrary airfoils allows greater design freedom 
in choosing incidence angle than exists with airfoils of a 
specified geometric family. Using a through-the-blade-row 
design approach and arbitrary airfoils, any reasonable distri- 


‘bution of incidence can be assumed and an apparently 


satisfactory solution obtained. Particularly with supersonic 
sections, the adequacy of the initial assumption should be 
tested by analyzing several sections in the cascade plane. 


The initial assumption made for. this rotor was a constant 
1.0 degree of incidence with respect to the suction surface. 
Because the relative inlet Mach number to the rotor is 


supersonic over most of the span, flow through the rotor. is 


controlled primarily by the wave pattern propagated upstream. 

A cascade~plane analysis using an inviscid, time-dependent 
caleulation technique shoved the initial rotor design to be 
about four percent deficient in the flow at both hub and tip. 
This corresponded to approximately one degree of inaidence 

with respect to the suction surface. The final design assumed» 
a constant 2.0 degrees of incidence with respect to the 

suction surface. The corresponding spanwise distribution of 
incidence with respect to the camber line is shown in Figure 11. 


The stator inlet Mach number varied from approximately sonic 


at the hub to 0.7 at the tip. Because of the sonic Mach number 


and high solidity at the hub conducive to choking and also the. 
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incidence tolerance presumably offered by leading edge sweep, 

a positive incidence of 5.0 degrees with respect to the 

camber line was assumed for the hub. Zero degrees of incidence | 
with respect to the camber line was assumed for the stator tip. 
This was selected on the basis of high cambered subsonic 
cascade sections favoring low to negative incidence angles and 
also the expectation of high incidence angles occurring here 

as the stage was throttled toward stall. The spanwise varia- 
tion from hub to tip was made linear and is also shown in 
Figure 11. No subsequent checks or adjustments were made to 
this distribution. 


10. DEVIATION ANGLE 


Deviation angle was predicted for both rotor and stator 
according to the method developed by NACA and described in 
Reference 8, Equations 269 and 271. A shape correction factor 
of 1.0 was used for both blade rows. However, the values 
predicted for the rotor were increased by 2.0 degrees at all 
radii, based upon examination of a few examples of potentially 
relevant recent data, The design deviation angle distribu- 
tions for both blade rows are presented in Figure 12. 


ll. AIRFOIL FILLET RADII 


The treatment of fillet radii at the rotor hub and both 
hub and tip of the stator merits mention because the approach 
was unconventional. For decades great attention has been 
paid to the design of fillets at the juncture cf an aircraft 
wing and fuselage. Poorly designed fillets, or no fillet at 
all, can lead to large interference drag penalties, Also, 
it is well known that in diffusing passages, boundary 
layer growth in corners is greater than on adjacent surfaces 
and boundary layer separation will generally occur first ina 
corner under an adverse pressure gradient. However, common 
practice in turbomachine design appears to favor the smallest 
possible fillets for aerodynamic reasons, mitigated only by 
manufacturing and structural desires for larger fillets. The 
fillet geometry chosen for this design is pattern:d after — 
the alreraft wing-fuselage model, although somewhat simplified. 
A radius of 0.25 inch was employed over most of the chord _ 
length; a rather large value in relation to most significant 
dimensions of this 17.0 inch diameter stage. The fillet 
radius declines smoothly from 0.25 inch to 0.06 inch over a 
distance of 6.78 inch approaching the leading and trailing edges 
of the.matar_and the leading edge of the stator, and over a | 
distance of 0.5 inch approaching the trailing edge of the stator. 
Note that the stator is hub-shrouded and the fillet treatment — 
applies to both platforms, — 


12. STRUCTURAL IMPACT ON AERODYNAMIC DESIGN 


A stress analysis of the steady state centrifugal and gas 
bending loads was performed for the ‘rotor as part of a mechanical 
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design and fabrication contract. Small adjustments to the 
section centroid locations with respect te the stacking axis 
were required and the leading and trailing edges approaching 
the hub platform were thickened before satisfactory stress 
levels were achieved. All of these adjustments were recycled 
through the aerodynamic design calculation so that the final 
aerodynamic design was fully consistent with the structural 
design. A peak combined stress level of about 68,000 psi 

was predicted to occur at the leading edge at a radius of 6.5 
inches at 20,000 revolutions per minute operating speed. 

The predicted rotor untwist distribution is shown in 

Figure 13. Correction for untwist was made by restaggering 
the cartesian airfoil manufacturing sections an equal and 
opposite amount. 


Vibratory tests of sample airfoils were also performed 
under the above mentioned contract in order to establish 
natural frequencies and mode shapes and to determine strain 
gage locations and dynamic stress operating limits. Three 
gage locations were chosen corresponding to points of 
maximum sensitivity for the three lowest order modes of 
vibration. These modes were 300 Hz (first bending) 1050 Hz 
(second bending), and 1400 Hz (first torsion). The first 
eight modes of vibration were mapped and stress ratios defined 
for the three gage locations. The eighth mode was a complex 
mode at 4900 Hz. An operating dynamic stress limit of 420,000 
psi was chosen. 
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SECTION IV 
FINAL DESIGN RESULTS 


1. STAGE DIMENSIONAL CHARACTERISTICS 


Stage Outer Diameter (constant) 17.00 inches , 
Rotor Inlet Hub/Tip Radius Ratio 0.312 
Number of Rotor Blades 20 
Number of Stator Blades 31 
Static Rotor Tip Clearance 0.027 
Design Point Rotor Running 
Clearance 0.016 
Rotor Aspect Ratio 1.320 
Stator Aspect Ratio 1.255 
2. DESIGN POINT SPECIFICATIONS 
Flow Rate 62.60 lb/sec 
Flow Per Unit Frontal Area 39.715 1b/sec/ft? 
Flow Per Unit Annulus Area 43.995 lb/sec/ft* 
Rotor Total Pressure Ratio 1.966 
Stage Total Pressure Ratio 1.912 
Rotor Tip Static Pressure Ratio 2.169 
Rotor Isentropic Efficiency 0.869 
Stage Isentropic Efficiency 0.830 
Inlet Corrected Tip Speed 1500 ft/sec 


3. AERODYNAMIC CHARACTERISTICS 


The details of the aerodynamic flow field throughout the 
stage are presented in the following pages of printout from 
the aerodynamic design computer program. The numbers and 
arrangement of computing stations are presented in Figure 14, 
For convenience, the computing stations defining blade row 
edges are also defined below. 


Station No. 11 
Station No. 16 
Station No. 18 
Station No. 23 


Roter leading edge 
Rotor trailing edge 
Stator leading edge 
Stator trailing edge 


The final streamwise distributions of total temperature through 


the rotor and radius-times~swirl-velocity through the stator 
are presented in Figures 15 and 16 respectively. The resultant 
streamwise static pressure distributions within both blade 

rows are shown in Figure 17. Figures 18 through 23 present a 
variety of spanwise parameter distributions of common interest. 
Rotor and stator relative inlet Mach number, diffusion factor, 
and total-pressure loss coefficient distributions are presented 
in Figures 18 through 20 respectively. Distributions of rotor- 
exit and stage-exit total pressure ratio, meridional velocity, 
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and isentropic efficiency are presented in Figures 21 through: 
23 respectively. The manufacturing dimensions of the annulus 


are shown in Figure 24, 
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4. ROTOR GEOMETRY 


Rotor airfoils were defined on all streamsurfaces and 
stacked close to the centroids cf the manufacturing sections 
on cartesian planes. The manufacturing sections were 
determined by mathematically passing a spline through common 
points on all streamsurface sections and determining the 
intersections of these splines with cartesian planes normal to 
the stacking axis. The computer program used to accomplish 
this was that described in Reference 7, previously cited. 
However, the interpolation technique and calculation of mechan- 
ical properties were taken from Reference 9 which discusses 
this subject more fully. The computer printout on the following 
pages presents detailed data on all streamsurface and manufac- 
turing sections. Except for the normalized data defining the 
streamsurface blade sections, all dimensions are in inches. 
On the first few pages appear sundry constants and a definition 
of the 21 streamsurfaces. The streamsurfaces are defined at 
eight axial locations which coincide with eight of the computing 
stations used for the aerodynamic design calculations. The 
origin for the axial locations of the stations is the same 
as was used for the aerodynamic analyses. The input data 
printout is completed with a table defining the geometry of 
each sectinn. Next are shown details of the 21 streamsurface 
sections. Only the "normalized" data have been reproduced; 
the equivalent dimensional data would be derived by scaling 
the nondimensional quantities by the meridional chord of the 
section. Finally, details of 13 manufacturing sections 
through the blade are shown. These plane sections, perpendicu~ 
lar to the stack axis, are _spaced 0,50 inch apart, and extend. 
slightly beyond the blade in both directions. The "ze 
coordinate is measured along the stack axis from the machine 
axis. The origin for the section coordinates is the stack | 
axis. The "X" direction is parallel to the machine axis, 
and "X" increases in the direction of flow. The "Y" direction 
is perpendicular to the "X" direction, and the "Y" coordinate 
decreases in the direction of rotation, "XS" and "YS" define 
the suction surface of the section, and "XP" and "YP" define 
‘the pressure surface, "XSEMI" and YSEMI" define the leading 
edge radius. The tralling edge is a straight line joining - 
the pressure and suction surfaces. The manufacturing 
coordinates shown in this report have not been corrected for 
blade untwist and represent the blade as it should exist when 
running at design speed. The actual machining was ‘accomplished 
with these coordinates rotated about the stacking axis an 
amount equal and opposite to the predicted angle of untwist 
as described in Section III.12. No corrections were made for 
any changes in canber-line shape. Figure 25 shows superimposed 
plots of developed streansurface sections. Figure 26 shows a 
similar view of the ,manufacturing sections. The larger change 
of section visible in Figure 26. is due to extrapolation of the 
airfoil beneath the hub ramp. 
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5. STATOR GEOMETRY 


The technique and ccnputer program used to define stator 
airfoils were the same as used for the rotor. The only 
significant differences were the stacking axis location and 
the number and spacing of manufacturing planes. The stator 
stacking axis was located near the trailing edge in order to 
minimize acute wall intersection angles in the aft portion 
of the passage. However, the stator twist turned out to be 
so slight that this is not a critical choice. Because of 
the reduced span (relative to the rotor) and the more 
complex leading edge shape, 11 manufacturing sections were 
employed spaced 0.375 inch apart. The computer printout 
on the following pages fully defines the stator airfoils 
and is identical in content and format to that shown for 
the rotor. Superimposed plots of the stacked streamsurface 
sections are shown in Figure 27 and of the manufacturing 
sections in Figure 28. Every other manufacturing section has 
been eliminated from the plot to improve the clarity of the 
figure. . 
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TOTAL PRESSURE LOSS PARAMETER 


TOTAL PRESSURE LOSS PARAMETER 
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FIGURE 1A. ROTOR L055 PARAMETER VS DIFFUSION FACTOR 
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FIGURE 18. STATOR LOSS PARAMETER VS DIFFUSION FACTOR 
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FRACTION OF MERIDIONAL CHORD FROM LEADING EDGE 


FIGURE 2. DESIGN DISTRIBUTIONS FOR INTRA=BLADE 
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FIGURE 10. STATOR SOLIDITY DISTRIBUTION 
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FIGURE 26, STACKED ROTOR CARTESIAN SECTIONS 
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APPENDIX A 


COMPUTER PROGRAM MODIFICATION 
FOR SUPERSONIC EXPANSION ANGLE 


One modification was incorporated into the computer 
program described in Reference 1. This modification allowed 
the user to specify the supersonic expansion angle 
directly, in lieu of specifying the relative flow angle at 
the shock. This option proved much easier to employ than the 
origianl options when neither the relative flow angles nor 
the cambers were known to within three or four degrees. 


The following card replacements are required: 


Change 1: Subprogram LOSS, Card LOSE 2223 
Replacement: 
25 FRDEL(L,J) = AA/RADIAN 


Change 2: Subprogram LOSS, Card LOSE 2224 
Replacement: 
FRDEL(L+1,J) = BB/RADIAN 


Change 3: Subprogram OUTPUT, Three continuation cards in 


format statement $1 (after JT. 2987) 
Replacement: 
X BSHPRESSURE 3X 16H DELTA B, IN - 3X 8H SOLIDITY 23X 
SH WHIRL 5X 
X16H DELTA B, IN = 3X 8HSOLIDITY/11X 7PROFILE 4X 
X16H LET TO SHOCK 33X 8HVELOCITY 3X 16H LET TO SHOCK//) 


This option is triggered by the integer 4 in columns 21-25 of 
Data Card number 5. Suitable descriptive corrections concern- 
ing tais option should be entered on pages D=3 (under Hesanne 
"Col. 2i-7i") and De10 Cunder “Card Type 20"). 
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